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( Dedicated to the 80th Anniversary of Steve Smale 

C"| ■ ABSTRACT Being neither commutative nor associative, Ein- 

stein velocity addition of relativistically admissible velocities gives 
rise to gyrations. Gyrations, in turn, measure the extent to which 
Einstein addition deviates from commutativity and from associa- 
tivity. Gyrations are geometric automorphisms abstracted from the 
relativistic mechanical effect known as Thomas precession. Thomas 
precession is a peculiar special relativistic space rotation that stems 
from Einstein velocity addition so that, in turn, Einstein velocity 
addition is algebraically and geometrically regulated by Thomas 
precession. Gyrations give rise to gyroalgebra and gyrogeometry, 
elegant nonassociative algebra and geometry in which groups and 
| vector spaces become gyrogroups and gyrovector spaces and Eu- 

clidean geometry becomes hyperbolic geometry. Indeed, gyrovec- 
tor spaces form the algebraic setting for the hyperbolic geometry 
. of Bolyai and Lobachevsky just as vector spaces form the algebraic 

setting for Euclidean geometry. As such, gyrations provide the 
missing link between Einstein's special theory of relativity and the 
hyperbolic geometry of Bolyai and Lobachevsky. 



In an exhaustive review of the vast literature on Thomas pre- 
cession, G.B. Malykin points out that many studies of Thomas 
precession are erroneous for various reasons. In particular, Ma- 
lykin's review reveals that there is no agreement among Thomas 
precession explorers on whether or not the Thomas precession an- 
gle of rotation and its generating angle have opposite signs. The 
mission of this article is to employ the study of gyrations and their 
geometry to the understanding of Thomas precession and, hence, 
to put Thomas precession misconceptions to rest. This article is 
dedicated to the 80th Anniversary of Steve Smale, whose work in 
geometric mechanics is well known. 
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1. Introduction 

It has been the lifetime desire of Steve Smale to improve our understanding of 
geometric mechanics [31 j . a theory that fits extraordinarily well into dynamical 
systems framework, as he explains in his 1967 survey article [30] . An overview of 
his involvement with geometric mechanics, without entering into technical details, 
is presented by Marsden in [21 . The impact of Einstein's work [25] led smale to 
remark in }3H p. 365] that relativity theory respects classical mechanics since "Ein- 
stein worked from a very deep understanding of Newtonian theory" . The present 
article on the hyperbolic geometric interpretation of the relativistic mechanical ef- 
fect known as Thomas precession is therefore dedicated to the 80th Anniversary of 
Steve Smale for his leadership and commitment to excellence in the field of geomet- 
ric mechanics. 

Thomas precession of Einstein's special theory of relativity is a physical realiza- 
tion of the abstract gyration. The latter, in turn, is an automorphism (defined in 
Def. |3]) that provides the missing link between Einstein's special theory of relativity 
and the hyperbolic geometry of Bolyai and Lobachevsky. Named after Llewellyn Hil- 
leth Thomas (1902-1992) who discovered its physical significance in 1926 [Tl 1331 IBIS] . 
Thomas precession is a special relativistic kinematic effect that regulates Einstein 
velocity addition both algebraically and geometrically [45] . In an exhaustive review 
of the vast literature on Thomas precession [20], G.B. Malykin emphasizes the im- 
portance of the frequency of the precession, pinpointing related erroneous results 
that are common in the literature. 

Accordingly, a study of Thomas precession in terms of its underlying hyperbolic 
geometry and elegant algebra is presented here in order to clarify the concept of 
Thomas precession. Thomas precession is an important special relativistic rotation 
that results from the nonassociativity of Einstein velocity addition and, hence, does 
not exist classically. Indeed, it was discovered in 1988 38 that Thomas precession 
regulates Einstein velocity addition, endowing it with a rich algebraic structure. 
As such, Thomas precession admits extension by abstraction, in which precession 
becomes gyration. The latter, in turn, gives rise to two new algebraic structures 
called a gyrogroup and a gyrovector space, thus introducing new realms to explore. 
The basic importance of gyrations is emphasized in gyrolanguage, where we prefix a 
gyro to any term that describes a concept in Euclidean geometry and in associative 
algebra to mean the analogous concept in hyperbolic geometry and nonassociative 
algebra. 

Gyrovector spaces turn out to form the algebraic setting for the hyperbolic ge- 
ometry of Bolyai and Lobachevsky just as vector spaces form the algebraic setting 
for Euclidean geometry [44j |45] [55] . This discovery resulted in the extension by ab- 
straction of Thomas precession into gyration, and in subsequent studies presented 
in several books [3U SSJ HH1 HS1 US EI] reviewed, for instance, in [55] and [M|. The 
hyperbolic geometric character of Thomas precession is emphasized here by the 
observation that it can be interpreted as the defect of a related hyperbolic triangle. 

Following Malykin's observations in [3D] , there is a need to demonstrate that our 
study of Thomas precession does lead to the correct Thomas precession angular 
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velocity. Accordingly, in this article we review the studies in [JS] 221 2H1 EH1 1101 IS] 
of both Thomas precession and its abstract version, gyration. Based on the review 
we derive the correct Thomas precession angular velocity, illustrated in Fig. 2] 

A boost is a Lorentz transformation without rotation |41j . The Thomas preces- 
sion angle e is generated by the application of two successive boosts with velocity 
parameters, say, u and v. The angle 9 between u and v is the generating angle of 
the resulting Thomas precession angle e, shown in Fig. |4] 

An important question about the Thomas precession angle e and its generating 
angle 9 is whether or not e and 9 have equal signs. According to Malykin [20], some 
explorers claim that e and 9 have equal signs while some other explorers claim that e 
and 9 have opposite signs. Malykin claims that these angles have equal signs while, 
in contrast, we demonstrate here convincingly that these angles have opposite signs. 
Our demonstration is convincing since it accompanies a focal identity, (jllOp . that 
interested explorers can test numerically in order to corroborate our claim that e 
and 9 have opposite signs. 

In order to pave the way to study Thomas precession we present the path from 
Einstein velocity addition to the gyroalgebra of gyrogroups and gyrations, and 
to the gyrogeometry that coincides with the hyperbolic geometry of Bolyai and 
Lobachevsky. We, then, demonstrate that the concept of Thomas precession in 
Einstein's special theory of relativity is a concrete realization of the abstract con- 
cept of gyration in gyroalgebra. 

A signed angle 9, —n < 9 < n, between two non-parallel vectors u and v in 
the Euclidean 3-space R 3 is positive (negative) if the angle 9 drawn from u to 
v is drawn counterclockwise (clockwise). The relationship between the Thomas 
precession signed angle of rotation, e, and its generating signed angle, 9, shown in 
Fig. |4l is important. Hence, finally, we pay special attention to the relationship 
between the Thomas precession signed angle of rotation and its generating signed 
angle, demonstrating that these have opposite signs. 



2. Einstein Velocity Addition and Scalar Multiplication 

Let (R 3 , +, •) be the Euclidean 3-space with its common vector addition, +, and 
inner product, •, and let 

(1) R 3 = {v e R 3 : ||v|| < c} 

be the c-ball of all relativistically admissible velocities of material particles, where 
c is the vacuum speed of light. 

Einstein velocity addition is a binary operation, ©, in the c-ball R 3 of all rel- 
ativistically admissible velocities, given by the equation, [45], [28] Eq. 2.9.2].[23] 
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p. 55], 0, 



(2) 



If 1 1 7u , v 
u©v = <^ u H v + -=• — '- (u-v u 



c 2 



„2 



1 + }±X { c 2 1 + 7, 



(u x (u x v)) 



for all u, v € R 3 , where 7 U is the gamma factor given by the equation 
(3) 7 V = ' 



cr 

Here u-v and ||v|| are the inner product and the norm in the ball, which the ball R 3 
inherits from its space R 3 , ||v|| 2 = vv = v 2 . Recalling that a nonempty set with 
a binary operation is called a groupoid, the Einstein groupoid (R 3 ,©) is called an 
Einstein gyrogroup. A formal definition of the abstract gyrogroup will be presented 
in Def. |H 

Einstein addition admits scalar multiplication ®, giving rise to the Einstein gy- 
rovector space (R 3 ,©,<g)). Remarkably, the resulting Einstein gyrovector spaces 
(R™,©,®) form the setting for the Cartesian-Beltrami-Klein ball model of hyper- 
bolic geometry, just as vector spaces form the setting for the standard Cartesian 
model of Euclidean geometry, as we will see in the sequel. 

Let fev be the Einstein addition of k copies of v € R™, that is fc<g5v = v©v . . . ©v 
(k terms). Then, 



(4) /c©v 



1 + MV_ (l _M 



The definition of scalar multiplication in an Einstein gyrovector space requires 
analytically continuing k off the positive integers, thus obtaining the following def- 
inition: 

Definition 1. (Einstein Scalar Multiplication; Einstein Gyrovector Spaces). 

An Einstein gyrovector space (R™, ffi, <S>) is an Einstein gyrogroup (R™, ffi) with scalar 
multiplication <S> given by 

<5) * V = ° 1+ Hiy + ( 1 - IMIV li^il = stanh( '' tani " ^Im, 

where r is any real number, r € R, v G R™, v / 0, and r©0 = 0, and with which 
we use the notation v<S>r = r®v. 
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In the Newtonian limit of large c, c — > oo , the ball R 3 expands to the whole of its 
space R 3 , as we see from ([1]), and Einstein addition © in R 3 reduces to the ordinary 
vector addition + in R 3 , as we see from @ and ©. 

Einstein addition @ of relativistically admissible velocities was introduced by 
Einstein in his 1905 paper [7] [8j p. 141] that founded the special theory of relativity. 
One has to remember here that the Euclidean 3- vector algebra was not so widely 
known in 1905 and, consequently, was not used by Einstein. Einstein calculated in 
[7] the behavior of the velocity components parallel and orthogonal to the relative 
velocity between inertial systems, which is as close as one can get without vectors 
to the vectorial version ([5} of Einstein addition. 

We naturally use the abbreviation uGv = uffi(- v) for Einstein subtraction, so 
that, for instance, v0v = 0, 0v = OGv = — v and, in particular, 

(6) e(uev) = euev 

and 

(7) Qu®(u®v) = v 

for all u, v in the ball R 3 , in full analogy with vector addition and subtraction in 
R 3 . Identity ^ is known as the automorphic inverse property, and Identity ([7]) 
is known as the left cancellation law of Einstein addition |48j . We may note that 
Einstein addition does not obey the naive right counterpart of the left cancellation 
law (0 since, in general, 

(8) (uffiv)Bv ^ u 

The seemingly lack of a right cancellation law for Einstein addition is repaired in 
(|77p by the introduction of a second binary operation, called Einstein coaddition, 
(|75p . which captures important analogies. 

Einstein addition and the gamma factor are related by the gamma identity, 

(9a) 7uev = 7 u 7v (l + ^r) 

which can be written, equivalently, as 

( 9b ) 7 eu ®v = 7 u 7v (i - -^r) 

for all u, v G R 3 . Here, (|9b[) is obtained from (|9"aj) by replacing u by 0u = u. 
A frequently used identity that follows immediately from (|3]) is 

v 2 llvll 2 7 2 -l 

(io) ^ = = 2^ 

c z cr 7^ 
and, similarly, useful identities that follow immediately from (0) are 

(Ha) u^ = _ 1+ 7^ 

c 7 u 7v 

and 

(lib) u ; v = 1 _7eu®v 

c 7u7v 
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implying 

( llc ) 7uev - 7u7v = -7 Q uev + 7u7v 

In Identity (jllcp the left-hand side seems to be more elegant, in form, than the 
right-hand side. Geometrically, however, the right-hand side of this identity is ad- 
vantageous over its left-hand side because the gamma factor 7q U V that appears on 
the right-hand side possesses a geometric interpretation. It is a geometric interpre- 
tation in hyperbolic triangles, called gyrotriangles, as explained in [5 1 1 Sec. 2.10]. 
To be more specific, we recall in Sec. 0] relevant results from [51] , 

Einstein addition is noncommutative. Indeed, in general, 

(12) uffiv ^ v®u 

u, v G M-c- Moreover, Einstein addition is also nonassociative since, in general, 

(13) (u©v)©w ^ uffi(vffiw) 
u, v, w <E Rg. 

It seems that following the breakdown of commutativity and associativity in 
Einstein addition some mathematical regularity has been lost in the transition from 
Newton's velocity vector addition in R 3 to Einstein's velocity addition @ in R 3 . 
This is, however, not the case since gyrations come to the rescue, as we will see 
in Sec. [5l Owing to the presence of gyrations, the Einstein groupoid (R 3 ,,©) has 
a grouplike structure [42j that we naturally call an Einstein gyrogroup [45J. The 
formal definition of the resulting abstract gyrogroup will be presented in Def. [H 
Sec. El 



3. Linking Einstein Addition to Hyperbolic Geometry 

The Einstein gyrodistance function, d(u, v) in an Einstein gyrovector space 
(R™, ©, <g>) is given by the equation 

(14) d(u,v) = ||uev|| 

u, v £ R™. We call it a gyrodistance function in order to emphasize the analogies it 
shares with its Euclidean counterpart, the distance function ||u — v|| in R™. Among 
these analogies is the gyrotriangle inequality according to which 

(15) ||uffiv|| < ||u||$||v|| 

for all u, v 6 R™. For this and other analogies that distance and gyrodistance 
functions share see [JSJ 28] . 

In a two dimensional Einstein gyrovector space (R^,©,®) the squared gyrodis- 
tance between a point x £ R 2 and an infinitesimally nearby point x + dx £ R 2 , 
dx = (dxi, dx 2 ), is defined by the equation [J51 Sec. 7.5] [321 Sec. 7.5] 

ds 2 = ||(x + dx)©x|| 2 

(16) 

= Edx\ + 2Fdx\dx 2 + Gdx\ + ... 
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where, if we use the notation r 2 = x\ + x 2 ,, we have 

E = c 2 



(c 2 — r 2 ) 2 

2 ^1^2 

(17) = c ^ 2 _ r 2)2 



G = c 2 



(c 2 



r2A2 



The triple (gn, <?i2, 922) = (-E 1 , F, G) along with 321 = 312 is known in differential 
geometry as the metric tensor gij [17) . It turns out to be the metric tensor of the 
Beltrami-Klein disc model of hyperbolic geometry [22] p. 220]. Hence, ds 2 in (fT^|) 
(fTT)) is the Riemannian line element of the Beltrami-Klein disc model of hyperbolic 
geometry, linked to Einstein velocity addition ([5J and to Einstein gyrodistance 
function [37] . 

The link between Einstein gyrovector spaces and the Beltrami-Klein ball model 
of hyperbolic geometry, already noted by Fock [9] p. 39], has thus been established 
in (|14[) - (|17l) in two dimensions. The extension of the link to higher dimensions is 
presented in 05J Sec. 9, Chap. 3], ggj Sec. 7.5] [46] Sec. 7.5] and [47]. For a brief 
account of the history of linking Einstein's velocity addition law with hyperbolic 
geometry see [2H1 p. 943]. 



4. Gyrotriangle, the Hyperbolic Triangle 

In this inspirational section we present recollections from [51] that intend to mo- 
tivate Thomas precession explorers to study the gyrostructure of Einstein addition 
and its underlying hyperbolic geometry. 

To be specific about the advantage in geometric interpretation of the use of 
Teuev over the use of 7 uffiv that appear in Identity (jllcj) we recall from [51] that 
Einstein addition admits scalar multiplication, eg), turning the Einstein gyrogroup 
(Rj?,®) into a corresponding Einstein gyrovector space as we will see in 

Seed 

Let U, V and W be the three vertices of a gyrotriangle UVW in an Einstein 
gyrovector space (R^,©,®), shown in Fig. [TJ Then, in full analogy with Euclidean 
geometry, the three sides of the gyrotriangle form the three gyrovectors 

u = ew®v 

(18) v = QW®U 

w = eu®v 

and the corresponding three side-gyrolengths of the gyrotriangle are 

u= ||u|| = ||eww|| 

(19) v = ||v|| = Hewwii 

w = ll w|| = \\ou®v\\ = ||eu®v|| 
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Figure 1. The gyrotriangle UVW in an Einstein gyrovector 
space (R™,®,<g>) is shown for n = 2. Its sides are presented 
graphically as gyrosegments that join the vertices. They form 
the gyrovectors u, v,w, side-gyrolengths, u, v, w, and gyroangles, 
a,j8,7< The gyrotriangle gyroangle sum is less than n, the differ- 
ence, (5 = 7r — (a + /3 + 7), being the gyrotriangular defect. 



and the side gamma factors of the gyrotriangle are, accordingly, 7 U , 7 V and 
(20) 7w = 7 Q uev 

Hence, by (HO]), the gamma factors 7 U , 7 V and 7q U V in (|llcl) can be interpreted 
geometrically as the gamma factors of the three sides of a gyrotriangle UVW. 

For later reference we recall here that the gyrotriangular defect 6 of the gyrotri- 
angle UVW is given in terms of the gyrotriangle side gamma factors by the equation 
[501 Theorem 2.32] [5TJ Theorem 6.11] 



(21) tan| = 



8 = V 1 + 2 7u7v7w ~ 7u ~ 7v ~ 7v 
1 + 7 U + 7v + 7w 



where 7w = 7eu©v 

It is the gamma identity (j9a]) that signaled the emergence of hyperbolic geometry 
in special relativity when it was first studied by Sommerfeld [32 and Varicak [52" l l53 | . 
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Historically, it formed the first link between special relativity and the hyperbolic 
geometry of Bolyai and Lobachevsky, recently leading to the novel trigonometry in 
hyperbolic geometry that became known as gyrotrigonometry, studied in [451 1461 

ill ill [501 EE]. 



5. The Gyrostructure of Einstein Addition 
Vector addition, +, in R 3 is both commutative and associative, satisfying 

u + v = v + u Commutative Law 

u + (v + w) = (u + v) + w Associative Law 

(22) 

for all u, v, w £ R 3 . In contrast, Einstein addition, ©, in R 3 is neither commutative 
nor associative. 

In order to measure the extent to which Einstein addition deviates from associa- 
tivity we introduce gyrations, which are maps that are trivial in the special cases 
when the application of © is associative. For any u,v £ R 3 the gyration gyr[u,v] 
is an automorphism of the Einstein groupoid (R 3 ,ffi) onto itself, given in terms of 
Einstein addition by the equation 

(23) gyr[u, v]w = ©(uffiv)ffi{uffi(vffiw)} 

for all u, v, w £ R 3 . 

We recall that an automorphism of a groupoid (S, ©) is a one-to-one map / of S 
onto itself that respects the binary operation, that is, f(a®b) = f(a)®f(b) for all 
a, 6 £ S. The set of all automorphisms of a groupoid (<!?, ©) forms a group, denoted 
Aut(S r , ©), where the group operation is given by automorphism composition. To 
emphasize that the gyrations of an Einstein gyrogroup (R 3 ,©) are automorphisms 
of the gyrogroup, gyrations are also called gyroautomorphisms. 

A gyration gyr[u, v], u, v £ R 3 , is trivial if gyr[u, v]w = w for all w £ R 3 . 
Thus, for instance, the gyrations gyr[0,v], gyr[v,v] and gyr[v, ©v] are trivial for 
all v £ R 3 , as we see from (J53]) and ©. 

Einstein gyrations, which possess their own rich structure, measure the extent to 
which Einstein addition deviates from both commutativity and associativity as we 
see from the gyrocommutative and the gyroassociative laws of Einstein addition in 
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the following list of identities [35113511351 

u©v = gyr[u, v](vffiu) 
u©(v®w) = (u©v)©gyr[u,v]w 
(ufflv)ffiw = u©(v©gyr[v,u]w) 
gyr[u©v,v] = gyr[u, v] 
gyr[u,v®u] = gyr[u,v] 
gyr[Gu, Gv] = gyr[u, v] 
(gyr[u, v])" 1 = gyr[v,u] 

(24) 

for all u, v, w € M%. 

Einstein addition is thus regulated by gyrations to which it gives rise owing 
to its nonassociativity, so that Einstein addition and its gyrations are inextricably 
linked. The resulting gyrocommutative gyrogroup structure of Einstein addition was 
discovered in 1988 [35J. Interestingly, gyrations are the mathematical abstraction 
of the relativistic mechanical effect known as Thomas precession [48l Sec. 10.3], as 
we will see in Sec. [9J 



1351 [5D1 [5T], each of which has a name: 

Gyrocommutative Law 
Left Gyroassociative Law 
Right Gyroassociative Law 
Gyration Left Loop Property 
Gyration Right Loop Property 
Gyration Even Property 
Gyration Inversion Law 



6. Gyrations 

Owing to its nonassociativity, Einstein addition gives rise in ([2"3"f to nontrivial 
gyrations 

(25) gyr[u,v] :R 3 C -> R 3 C 

for any u, v G M%.. Gyrations, in turn, regulate Einstein addition, endowing it with 
the rich structure of a gyrocommutative gyrogroup that will be formalize in Sec. [7] 

The gyration equation is expressed in (|23p in terms of Einstein addition. Ex- 
pressing it explicitly, in terms of vector addition and vector scalar product rather 
than Einstein addition, we obtain the equation 

r i Au + Bv 

(26) gyr[u, vjw = w + 



where 



D 



1 7 2 1 

A = XTT ^v ~ l)(u-w) + -7u7 v (v-w) 

2 7 2 7 2 



(27) 



B = — 



c 4 (7u + l)(7v + l) 
1 7- 



c 2 7 V + 



j{7u(7v + l)(u-w) + ( 7u - l) 7v (v-w)} 



D = 7u7v(l + ^j?) + 1 = 7uev + 1 > 1 
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for all u, v,w € R 3 . Clearly, the domain of u and v in (|27|l must be restricted to 
the open ball R 3 to insure the reality of the Lorentz factors 7 U and j v . In contrast, 
however, the domain of w need not be restricted to the ball. 

Allowing w e R 3 D R 3 in ([5BJ -([2"7)l. that is, extending the domain of w from R 3 
to R 3 , gyrations gyr[u, v] are expendable from self-maps of R 3 to linear self-maps 
of R 3 for all u, v e R 3 . Indeed, 

(28) gyr[u, v](r 1 wi +r 2 w 2 ) = ^gyr^vjwi + r 2 gyr[u, v]w 2 

for all u, v e R 3 , w e R 3 and r 1 ,r 2 e R. 

In each of the three special cases when (i) u = 0, or (ii) v = 0, or (hi) u and v 
are parallel in R 3 , u||v, we have Au + Bv = 0, so that in these cases gyr[u, v] is 
trivial. Thus, we have 

gyr[0,v]w = w 

(29) gyr[u, 0]w = w 

gyr[u, v]w — w, u||v 

for all u, v e R 3 C R 3 and all w e R 3 . 
It follows from flU - (J2TJ) that 

(30) gyr[v,u](gyr[u,v]w) = w 

for all u, v 6 R 3 , w 6 R 3 , so that gyrations are invertible linear maps of R 3 , the 
inverse, gyr _1 [u, v], (j24|) . of gyr[u, v] being gyr[v, u]. We thus obtain from (f30|) the 
gyration inversion property in (|24|) , 

(31) gyr _1 [u,v] = gyr[v,u] 
for all u, v e R 3 . 

Gyrations keep the inner product of elements of the ball R 3 invariant, that is, 

(32) gyr[u, v]a-gyr[u, v]b = ab 

for all a, b, u,v G R 3 . Hence, in particular, gyr[u, v] is an isometry of R 3 , keeping 
the norm of elements of the ball R 3 invariant, 

(33) ||gyr[u,v]w|| = ||w|| 

Accordingly, for any u,v € R 3 , gyr[u,v] represents a rotation of the ball R 3 about 
its origin. 

The invertible self-map gyr[u,v] of R 3 respects Einstein addition in R 3 , 

(34) gyr[u, v] (a©b) = gyr[u, v]a©gyr[u, v]b 

for all a, b, u,v g R 3 , so that, by (|3T|) and , gyr[u,v] is an automorphism of 
the Einstein groupoid (R 3 ,©). 

7. Gyrogroups 

Taking the key features of the Einstein groupoid (R 3 ,ffi) as axioms, and guided 
by analogies with groups, we are led to the formal gyrogroup definition in which 
gyrogroups turn out to form a most natural generalization of groups. Definitions 
related to groups and gyrogroups thus follow. 
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Definition 2. (Groups). A groupoid (G, +) is a group if its binary operation 
satisfies the following axioms. In G there is at least one element, 0, called a left 
identity, satisfying 

(Gl) 0+a=a 

for all a £ G. There is an element G G satisfying axiom (Gl) such that for each 
a G G there is an element —a € G, called a left inverse of a, satisfying 

(G2) -a + a = 

Moreover, the binary operation obeys the associative law 
(G3) (a + b)+c = a+(b + c) 

for all a, b, c G G. 

Groups are classified into commutative and noncommutative groups. 

Definition 3. (Commutative Groups). A group (G, +) is commutative if its 
binary operation obeys the commutative law 

(G6) a + b = b + a 

for all a, b G G. 

Definition 4. (Gyrogroups). A groupoid (G, ffi) is a gyrogroup if its binary 
operation satisfies the following axioms. In G there is at least one element, 0, called 
a left identity, satisfying 

(Gl) 0ffia = a 

for all a G G. There is an element G G satisfying axiom (Gl) such that for each 
a G G there is an element Qa G G, called a left inverse of a, satisfying 

(G2) Qa®a = . 

Moreover, for any a, b, c G G there exists a unique element gyr[a, b]c G G such that 
the binary operation obeys the left gyroassociative law 

(G3) a®(b®c) = (a®b)®gyr[a,b]c. 

The map gyr[a, b] : G — > G given by c ^ gyr[a, &]c is an automorphism of the 
groupoid (G, ®), that is, 

(G4) gyr[a, b] G Aut(G, ©) , 

and the automorphism gyr[a, b] of G is called the gyroautomorphism, or the gyration, 
of G generated by a,b G G. The operator gyr : G x G — > Aut(G, ©) is called the 
gyrator of G. Finally, the gyroautomorphism gyr[a, b] generated by any a, b G G 
possesses the left loop property 

(G5) gyr[a, b] = gyr[a©&, b] . 

The gyrogroup axioms (G1)-(G5) in Def. |4]are classified into three classes: 

(1) The first pair of axioms, (Gl) and (G2) in Def. [H is a reminiscent of the 
group axioms (Gl) and (G2) in Def. [2j 

(2) The last pair of axioms, (G4) and (G5) in Def. 01 presents the gyrator 
axioms. 
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(3) The middle axiom, (G3) in Def. 01 is a hybrid axiom linking the two pairs 
of axioms in items (1) and (2). 

As in group theory, we use the notation aQb = a©(0&) in gyrogroup theory as 
well. 

In full analogy with groups, gyrogroups are classified into gyrocommutative and 
non-gyrocommutative gyrogroups. 

Definition 5. (Gyrocommutative Gyrogroups). A gyrogroup (G, ©) is gyro- 
commutative if its binary operation obeys the gyrocommutative law 

(G6) a®b = gyr[a,6](6©a) 

for all a, b G G. 

While it is clear how to define a right identity and a right inverse in a gyrogroup, 
the existence of such elements is not presumed. Indeed, the existence of a unique 
identity and a unique inverse, both left and right, is a consequence of the gyrogroup 
axioms, as the following theorem shows, along with other immediate results. 

Theorem 6. (First Gyrogroup Properties). Let (G, ©) be a gyrogroup. For 
any elements a,b,c,x € G we have the following results: 

(1) If a©6 = a©c, then b = c (general left cancellation law; see item (9) below). 

(2) gyr[0, a] = I for any left identity in G. 

(3) gyr[cc, a] = I for any left inverse a; of a in G. 

(4) gyr[a, a] = I 

(5) There is a left identity which is a right identity. 

(6) There is only one left identity. 

(7) Every left inverse is a right inverse. 

(8) There is only one left inverse, ©a, of a, and ©(©a) = a. 

(9) The Left Cancellation Law: 

(35) ©affi(affife) = b 

(10) The Gyrator Identity: 

(36) gyr[a, b]x = e(a®&)©{a®(6®x)} 

(11) gyr[a,6]0 = 0. 

(12) gyr[a,6](ea;) = Qgyi[a,b]x. 

(13) gyr[a,0]=J. 

Proof. 

(1) Let x be a left inverse of a corresponding to a left identity, 0, in G. We 
have x(B(a(Bb) = x®(affic), implying (x®a)(Bgyv[x, a]b — (x(Ba)®gyr{x, a]c 
by left gyroassociativity. Since is a left identity, gyr[a;,a]6 = gyr[a;,o]c. 
Since automorphisms are bijective, b = c. 

(2) By left gyroassociativity we have for any left identity of G, a@x = 
0(B(a(Bx) = (0©a)©gyr[0, a]x = a©gyr[0, a]x. Hence, by item [1] above we 
have x — gyr[0, a]x for all x £ G so that gyr[0, a] = I. 
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(3) By the left loop property and by item[2]above we have gyr[a;, a] = gyr[x®a, a] = 
gyr[0, a] = I. 

(4) Follows from an application of the left loop property and item [5] above. 

(5) Let x be a left inverse of a corresponding to a left identity, 0, of G. Then 
by left gyroassociativity and item|3]above, x©(a©0) = (x©a)©gyr[a;, a]0 = 
OffiO = = x@a. Hence, by (1), a®0 = a for all a e G so that is a right 
identity. 

(6) Suppose and 0* are two left identities, one of which, say 0, is also a right 
identity. Then = 0*©0 = 0*. 

(7) Let x be a left inverse of a. Then x®{a(&x) = (x©a)©gyr[x, a]x = 0(Bx = 
x — x(B0, by left gyroassociativity, (G2) of Def. 2] and items [31 El El above. 
By item [T] we have a(Bx = so that a; is a right inverse of a. 

(8) Suppose x and y are left inverses of a. By item [7] above, they are also right 
inverses, so a®x = = affiy. By item[TJ x = y. Let Go be the resulting 
unique inverse of a. Then 0«©o = so that the inverse 0(0<z) of ©a is a. 

(9) By left gyroassociativity and by [3] we have 

(37) ©offi(affio) = (ea©a)©gyr[Ga, a]b = b 

(10) By an application of the left cancellation law in item [!|] to the left gyroas- 
sociative law (G3) in Def. 2] we obtain the result in item [TU] 

(11) We obtain item [HI from [TQl with x = 0. 

(12) Since gyr[a, b] is an automorphism of (G, ©) we have from [TT1 

(38) gyr[a, b] (Gx)©gyr[a, b]x = gyr[a, b] (Qx®x) = gyr[a, b}0 = 

and hence the result. 

(13) We obtain item [13] from [10] with 6 = 0, and a left cancellation, item [9] 

Einstein addition admits scalar multiplication, giving rise to gyrovector spaces. 
Studies of gyrogroup theory and gyrovector space theory along with applications in 
hyperbolic geometry and in Einstein's special theory of relativity are presented in 

^ISSlllSlllSlEBEIl- 

We thus see in this section that Einsteinian velocity addition, ©, in of relativis- 
tically admissible velocities gives rise to the gyrocommutative gyrogroup (]Rj?,ffi), 
just as Newtonian velocity addition, +, in R 3 of classical, Newtonian velocities 
gives rise to the commutative group (R 3 , +). Newtonian velocity addition, in turn, 
is given by the common vector addition in R 3 . Clearly, it is owing to the presence of 
nontrivial gyrations that gyrogroups are generalized groups. Accordingly, gyrations 
provide the missing link between the common vector addition and Einstein velocity 
addition. 

8. The Euclidean and Hyperbolic Lines 

As shown in Figs. [5]-[31 we introduce Cartesian coordinates into R ra in the usual 
way in order to specify uniquely each point P of the Euclidean n-space R™ by an 
n-tuple of real numbers, called the coordinates, or components, of P. Cartesian 
coordinates provide a method of indicating the position of points and rendering 
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m AB =A + {-A + B)\ 
d(A, B) = \\A-B\\ 
-•1* d(A i m A B ) = d(B,m A B ] 

d(A, P) + d{P, B) = d{A, B) 



A + {-A + B)t 



— oo < t < oo 



-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.6 



Figure 2. The Cartesian coordinates for the Euclidean plane 
R 2 , (xi,x 2 ), x\ + x\ < oo, are shown. The points A and B in 
the Euclidean plane R 2 are given, with respect to these Cartesian 
coordinates, by A = (-0.60,-0.15) and B = (0.18,0.80). The 
distance function, d(A, B) = \\A— B\\ , and the equation of the line 
through A and B are shown along with the triangle equality for 
a generic point P on the line. The midpoint m A B of the segment 
AB is reached when the line parameter is t = 1/2. 



graphs on a two-dimensional Euclidean plane R 2 and in a three-dimensional Eu- 
clidean space R 3 . 

As an example, Fig. [2] presents a Euclidean plane R 2 equipped with a Cartesian 
coordinate system E. The position of points A and B and their midpoint Mab 
with respect to E are shown. 

The set of all points 
(39) A + (-A + B)t 

t € M, forms a Euclidean line. The segment of this line, corresponding to 1 < t < 1, 
and a generic point P on the segment, are shown in Fig. [2] Being collinear, the 
points A,P and B obey the triangle equality d(A,P) + d(P,B) = d(A,B), where 
d(A, B) = || — A + B\\ is the Euclidean distance function in R". 

Fig. [2] demonstrates the use of the standard Cartesian model of Euclidean ge- 
ometry for graphical presentations. In a fully analogous way, Fig. [3] demonstrates 
the use of the Cartesian-Bcltrami-Klein model of hyperbolic geometry; see also [51] 
Figs. 2.3-2.4]. 

Now, let A,Bd M™ be two distinct points of the Einstein gyrovector space 
(R™,©,(g)), and let t G R be a real parameter. Then, in full analogy with the 
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Figure 3. The Cartesian coordinates for the unit disc in the 
Euclidean plane R 2 , (xi, £ 2 ), x\ + x\ < 1, arc shown. These, inside 
the disc, are mow considered the Cartesian coordinates for the 
Einstein gyrovector plane (R 2 ,©,<gi). The points A and B in the 
Einstein gyrovector plane (M 2 , ©, ®) are given, with respect to these 
Cartesian coordinates, by A = (-0.60, -0.15) and B = (0.18, 0.80). 
The gyrodistance function, d(A, B) = \\A<3B\\, and the equation of 
the gyroline through A and B are shown along with the gyrotriangle 
equality for a generic point P on the gyroline. The gyromidpoint 
m A B of the gyrosegment AB, reached when the gyroline parameter 
is t = 1/2, shares obvious analogies with its Euclidean counterpart, 
the midpoint in Fig. [2] 

Euclidean line (|39|) . the graph of the set of all points, Fig. 02 



f e 1, in the Einstein gyrovector space (K",©, ®) is a chord of the ball R™. As 
such, it is a geodesic line of the Cartesian-Beltrami-Klein ball model of hyperbolic 
geometry. The geodesic line (|40|) is the unique geodesic passing through the points A 
and B. It passes through the point A when t = and, owing to the left cancellation 
law, (f35|) . it passes through the point B when t=l. Furthermore, it passes through 
the midpoint Ma,b of A and B when t = 1/2. Accordingly, the gyrosegment that 
joins the points A and B in Fig. [3] is obtained from gyroline (HPT) with < t < 1. 



It is owing to the gyrocommutative law, ([24]) . of Einstein addition that Thomas 
precession of Einstein's special theory of relativity is recognized as a concrete ex- 
ample of the abstract gyrogroup gyration in Def. |4l Accordingly, the gyrogroup 



(40) 



A®(oA®B)®t 



9. Thomas Precession 
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gyration is an extension by abstraction of the relativistic mechanical effect known 
as Thomas precession. 

The gyrocommutativc law of Einstein velocity addition was already known to 
Silberstein in 1914 [29 in the following sense: According to his 1914 book, Sil- 
berstein knew that the Thomas precession generated by u, v £ Rj? is the unique 
rotation that takes v©u into u©v about an axis perpendicular to the plane of u 
and v through an angle < tt in M 3 , thus giving rise to the gyrocommutative law. 
However, obviously, Silberstein did not use the terms "Thomas precession" and "gy- 
rocommutative law" . These terms have been coined later, respectively, (i) following 
Thomas' 1926 paper [35], and (ii) in 1991 [HIS], following the discovery of the 
accompanying gyroassociative law of Einstein addition in 1988 38 . 

A description of the 3-space rotation, which since 1926 is named after Thomas, 
is found in Silberstein's 1914 book In 1914 Thomas precession did not have 
a name, and Silberstein called it in his 1914 book a "certain space-rotation" [29l 
p. 169]. An early study of Thomas precession, made by the famous mathematician 
Emile Borel in 1913, is described in his 1914 book [5] and, more recently, in [33] . 
According to Belloni and Reina Q] , Sommerfeld's route to Thomas precession dates 
back to 1909. However, prior to Thomas discovery the relativistic peculiar 3-space 
rotation had a most uncertain physical status (54] p. 119]. The only knowledge 
Thomas had in 1925 about the peculiar relativistic gyroscopic precession [15], how- 
ever, came from De Sitter's formula describing the relativistic corrections for the 
motion of the moon, found in Eddington's book [5], which was just published at 
that time g3J Sec. 1, Chap. 1]. 

The physical significance of the peculiar rotation in special relativity emerged in 
1925 when Thomas relativistically re-computed the precessional frequency of the 
doublet separation in the fine structure of the atom, and thus rectified a missing 
factor of 1/2. This correction has come to be known as the Thomas half @], pre- 
sented in ([95]) . Thomas' discovery of the relativistic precession of the electron spin 
on Christmas 1925 thus led to the understanding of the significance of the relativis- 
tic effect which became known as Thomas precession. Llewellyn Hilleth Thomas 
died in Raleigh, NC, on April 20, 1992. A paper [3] dedicated to the centenary of 
the birth of Llewellyn H. Thomas (1902-1992) describes the Bloch gyrovector of 
quantum information and computation. 

Once recognized as gyration, it is clear that Thomas precession owes its exis- 
tence solely to the nonassociativity of Einstein addition of Einsteinian velocities. 
Accordingly, Thomas precession has no classical counterpart since the addition of 
classical, Newtonian velocities is associative. 

It is widely believed that special relativistic effects are negligible when the ve- 
locities involved are much less than the vacuum speed of light c. Yet, Thomas 
precession effect in the orbital motion of spinning electrons in atoms is clearly ob- 
served in resulting spectral lines despite the speed of electrons in atoms being small 
compared with the speed of light. One may, therefore, ask whether it is possible 
to furnish a classical background to Thomas precession |19) . Hence, it is important 
to realize that Thomas precession stems from the nonassociativity of Einsteinian 
velocity addition, so that it has no echo in Newtonian velocities. 
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In 1966, Ehlers, Rindler and Robinson [5] proposed a new formalism for dealing 
with the Lorentz group. Their formalism, however, did not find its way to the 
mainstream literature. Therefore, thirty three years later, two of them suggested 
considering the "notorious Thomas precession formula" (in their words |27[ p. 431]) 
as an indicator of the quality of a formalism for dealing with the Lorentz group. 
The idea of Rindler and Robinson to use the "notorious Thomas precession formula" 
as an indicator works fine in the analytic hyperbolic geometric viewpoint of spe- 
cial relativity, where the ugly duckling of special relativity, the "notorious Thomas 
precession formula" , becomes the beautiful swan, the "gyrostructure" , of special 
relativity and its underlying analytic hyperbolic geometry. The abstract Thomas 
precession, called gyration, is now recognized as the missing link between classical 
mechanics with its underlying Euclidean geometry and relativistic mechanics with 
its underlying hyperbolic geometry. 



10. Thomas Precession Matrix 



For any two vectors a, b £ IR 3 , a = (01,02,03), etc., determined by their com- 
ponents with respect to a given Cartesian coordinate system, we define the square 
3x3 matrix f2(a, b) by the equation 

(ai&i ai& 2 aib 3 \ (a\b\ a 2 b\ a 3 &i N 

a 2 bi a 2 b 2 a 2 b 3 + a x b 2 a 2 b 2 a 3 b 2 

a 3 h a 3 b 2 a 3 b 3 ) \a\b 3 a 2 b 3 a 3 b 3/ 

or, equivalently, 



(42) 


[ ° 

fi(a, b) = - -u 3 
\ u 2 


W3 



-ui2 



where 








(43) 


w = (cji, UJ 2 , L0 3 


,)=a 


x b 


Accordingly, 








(44) 


fi(a, b)x = (a x b) x x = 


= -a(b 


■x) + 


for any xell 


Hence, 






(1) 0(a,b) 

(2) and 


= if and only if a x b = 






(45) 


fi(a, b)(a x 


b) = 





(3) and, for Q = 0(a,b), 

(46) !l 3 = -(axb) 2 !! 

The matrix Q — Q(u, v) can be used to simplify the presentation of both Einstein 
addition u®v and its associated gyration gyr[u, v], 

(47) uev= ^_j u + v _^_Ji_o u 
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(48) gyr[u, v] = I + ail + fitf 

where / is the 3x3 identity matrix, and where 



(49) 





_^ /u Iv 

c 4 (l + 7u)(l + 7 v )(l + 7uev) 



1 



satisfying a < 0, j3 > 0, and 

(50) a 2 + [u 2 v 2 - (u-v) 2 ]/? 2 -2/3 = 

for all u, v e R 3 . 



The gyration matrix gyr[u, v] in (|48|) satisfies the cubic equation 



(51) 



gyr 3 [u, v] - trace(gyr[u, v] )gyr 2 [u, v] 



+ trace(gyr[u, v])gyr[u, v] — / = 



called the trace identity. 

The trace identity (|51[) characterizes 3x3 matrices that represent proper rotations 
of the Euclidean 3-space R 3 about its origin. 

The matrix representation of gyr[u, v] in R 3 relative to an orthonormal basis is 
thus an orthogonal 3x3 matrix with determinant 1. It follows from (|45l) and (14"5|) 
that 

(52) gyr[u,v](u x v) = u x v 

so that the vector u x v lies on the rotation axis of the gyration gyr[u, v]. 

Interesting studies of the trace identity, using analysis, algebra and geometry is 
found in an elementary form in |16) and in a more advanced form in [131 112[ 1111 110] . 



Let £", £' and £ be three inertial frames in the Euclidean 3-space R 3 with re- 
spective spatial coordinates (x", y"), (x' , y') and (x, y). The third spatial coordinate 
of each frame is omitted for simplicity. Accordingly, these are shown in Fig. [4] in R 2 
rather than R 3 . Frame £" moves with velocity v £ R 3 , without rotation, relative 
to frame £' which, in turn, moves with velocity u £ R 3 , without rotation, relative 
to frame £. The angle between u and v is 9, shown in Fig. 21 satisfying 



so that, by (fTO]) . 



Observers at rest relative to £ and observers at rest relative to £' agree that 
their coordinates (x, y) and (x',y') are parallel. Similarly, observers at rest relative 
to £' and observers at rest relative to £" agree that their coordinates (x',y') and 
(x" , y") are parallel, as shown in the left part of Fig. S) 



11. Thomas Precession Graphical Presentation 



(54) 



1 



— 7u7v u-v = 
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y y 




E S 



Figure 4. In the Euclidean plane R 2 an inertial frame E" moves 
uniformly, without rotation, with velocity v € R 2 relative to in- 
ertial frame E'. The latter, in turn, moves uniformly, without 
rotation, with velocity u £ R 2 relative to inertial frame E. Ow- 
ing to the presence of Thomas precession, the inertial frame E" 
moves uniformly, with rotation angle e, with a composite velocity 
relative to the inertial frame E. Is the composite velocity of E" 
relative to E uffiv or vffiu? The answer is: neither; see (|75|) . The 
Thomas precession signed angle e, — tt < e < 7T, turns out to be 
the unique rotation angle with rotation axis parallel to u x v in 
M.j? that takes vffiu into uffiv according to the gyrocommutative 
law uffiv = gyr[u, v](vffiu). Being related by (|57)l . the Thomas 
precession signed angle e and its generating signed angle 9 from u 
and v have opposite signs, illustrated graphically in Figs. ??-??. 



Counterintuitively, if 9 and 9 ^ tt, observers at rest relative to E and 
observers at rest relative to E" agree that their coordinates are not parallel. Rather, 
they find that their coordinates are oriented relative to each other by a Thomas 
precession angle e, < e < tt, as shown in the right part of Fig. 01 

Let u and v be two nonzero vectors in the ball Rj?. By the gyrocommutative law 
in ([24l . the gyration gyr[u, v] takes the composite velocity vffiu into uffiv. Indeed, 
gyr[u, v] is the unique rotation with rotation axis parallel to u x v that takes v©u 
into uffiv through the gyration angle e, < e < ir. We call e the Thomas precession 
(or, rotation) angle of the gyration gyr[u,v], and use the notation 

(55) e = Zgyr[u,v] 

Accordingly, the Thomas precession angle e = Zgyr[u, v] generated by u, v £ R^, 
shown in the right part of Fig. 01 satisfies the equations 

(uffiv) • (vffiu) 

cos e = n fi^ 

||uffiv|| 2 

(56) 

| (uffiv) x (vffiu) || 



sin e = ± - 
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Figure 5. A graphical presentation of the cosine of the Thomas 
precession angle e, cose, (|57|) , as a function of the angle 9 between 
its two generating relativistically admissible velocities u, v <E M 3 
for several values of k, k being a function, (|58jl. of 7 U and 7 V . 



The Herculean task of simplifying (|5o) was accomplished in [351 I3S1 ESI H2]) 
obtaining 

(fc + cos(9) 2 - sin 2 9 

c °s e = ^ o— 

(A; + cos6>) 2 -I- sin 6 1 

(57) 

-2(k + cos 6) sin 6* 

sm e = — 

(fc + cos(9) 2 +sin 2 6» 

where 9, < 9 < 2ir, is the angle between the vectors u,v £ M 3 , forming the 
horizontal axes in Figs.[5Hni and where k, k > 1, is a velocity parameter given by 
the equation 



(58) k z = 



2 _ 7u + 1 7v + 1 



7 U - 1 7v - 1 



The parameter k approaches 1 when both ||u|| and ||v|| approach c. We clearly 
have the limits 

lim cos e = cos 9 

(59) . 

lim sm e = — sm 9 

k^l 

for < 9 < 2vr, 9 ^ n, seen in Figs.[5]and[6] 
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TT 27T 

Figure 6. A graphical presentation of the negative sine of the 
Thomas precession angle e, — sine, (|57[) . as a function of the angle 
9 between its two generating relativistically admissible velocities 
u, v G IR 3 for several values of k, k being a function, (|58|) . of 7 U 
and 7 V . 



Figs. [5] and [6] present graphically cose and — sine as functions of 9 for several 
values of k. As expected, the graphs in these figures show that for all values of the 
parameter k, k > 1, Thomas precession angle e vanishes when 9 = 0, when 9 = n, 
and again, when 9 = 2tt. In the limit of high relativistic speeds approaching the 
vacuum speed of light c, ||u||, ||v|| — > c, the parameter k approaches unity, k — > 1, 
and e — > —9 for all 9 in the punctured interval [0, 7r)(J(7r, 2tt\. The punctured 
interval is the union of the two connected intervals [0, tt) and (tt, 2tt] which is the 
closed connected interval [0, 2tt] from which the point tt has been removed. Thus, 
there is no Thomas precession angle tt; see also ([70| . 

The extension by abstraction of Thomas precession into gyrations enables the de- 
velopment of techniques that explain the non-existence of a gyration whose rotation 
angle is tt; see the Gyration Exclusion Theorem in (48j Theorem 3.36]. 

As we see from Figs.[5]and[6l the variation of e for < 9 < 2ir is over the interval 
[0, 2tt] punctured by a fc-dependent subinterval centered at e = tt. 
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It is interesting to derive cos % and sin | from ([5 



e / 1 + cos e k + cos 9 

cos - = ± \ 



2 ^(fc + costf) 2 +sin 2 6» 

(60) 



e / 1 — cos e 



^ (k + cos#) 2 + sin 2 



Following Fig.01 the ambiguous signs in (|60|) are selected such that cos § > while 
sin | and sin# have opposite signs. 

12. Thomas Precession Angle 

Thomas precession gyr[u, v] in (J3HJ) can be recast into a form familiar as the 
representation of a rotation about an axis by an angle e, 



(61) gyr[u,v] 



J + sme-^ + (l-cose) ' , uj ^ 
. I, uj e =0 
where u, v <E R^, and where e is the Thomas precession angle shown in Fig.|4] 
Comparing (T5TT) with (14"5)) , we see that 

sine = a(u, x)ujq 

(62) 

1 — cos e = /3(u, v)wg 

and 

UJg = ±||u X v|| 

= N|||v||sin0 

(63) 



= c — sin & 

7u7v 

where the ambiguous sign is selected such that ojq and sin have equal signs. 

It follows from (j6"2"]) - (|63l) . and from the definition of a(u, v) and /3(u, v) in (14"!)]) 
that 



(64) 



co S e = l- (7 "- 1)( V 1) ^ 2 ^ 

7u©v + 1 



y^TT^l + (7u - l)(7v ~ 1) . a 

— 1 sm 



7uQ 

Following © - (HI]) and ((54]) we have 



(65a) 7 uev = 7 U 7 V + ~ W7v ~ 1 cos ' 

and 



(65b) 7euev = 7u7 v - VSS - Wl* ~ 1 cos( 
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so that, by flM) - flUHg 

cos, = 1 (7 u -lK7v-l^_ sin2 Q 

1 + 7 u 7v + yjll - V^v - !cos6> 

(66) 

(7„-l)(7v-l)(fe + cos0) . 
sm e = . . sin w 

1 + 7 u 7v + V72 - !V7v - lc °s0 
where k > 1 is given by (|58|). 

The special case when u and v have equal magnitudes is required for later refer- 
ence related to Fig. [4] In this special case 7 U = 7v , so that e in (|66[) reduces to e s 
given by 

_ ( 7v -l) 2 sin 2 g 
C0S£s l + 7 2 + (73-l)cos0 

(67) 

( 7 g-l) + ( 7v -l) 2 cosfl . 
SmC ' = - l + 7 2 + (7 2 -l)cos. 



Solving (1651) for cos 6 1 we obtain the equations 

7uev " 7u7v _ -7 Q uev + 7 u 7v 



cos 9 



(68) sin 2 = 1- cos 2 6 

= 1 - 7u - 7 2 ~ 7uev + 2 7u7v7uev 
(7S-1)(7 2 -1) 
The substitution of (|B5]l into ([Ml) gives 

1 

cos e = — — r 

(7u + l)(7v + l)(7u®v + l) 

(69) _ r ,2,2,2 

v ; x {-7u7v7uev + 7 U + 7v + 7 u ©v 

+ 7u7v + 7u7u©v + 7v7 u ®v + 7u + 7v + 7 u ©vl 
so that, finally, we obtain the elegant expression 

r<?n\ i _i_ (1 + 7 U +7 V +7 U ffiv) 2 . n 

(70) l + cose= — -j- — — - r>0 

(l + 7u)(l + 7v)(l + 7uev) 

which agrees with McFarlane's result, cited in [351 Eq. (2.10.7]. It implies that e ^ ir 
for all u, v G R;?; and that 



(71) cos - = / + C0Se = l + 7u+7v + 7uev 

2 V 2 v^v/i + T^v/ITt^ 1 + 7u©v 

Finally, we also have the elegant identity 

1 + 2 7 u 7v7u©v - 7u - 7v - 7ug 



(72) tan 2 



l + cose/ (1 + 7u + 7v + 7u©v) 2 
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Hence, (|55p . the Thomas precession angle e = Zgyr[u, v] in Fig. 2] is given by the 
equation 

/ 73 x tan 2 ^gyr[u, v] = 1 + 27u7v7uev ~ ll ~ 7v ~ 7u @ v 

2 (i +7u + 7v + 7 )2 



Noting the gyration even property in (|24|) and replacing u by Gu in (|73[) , we obtain 
the equations 

tan 2 Zgyr[u ' Qv] = tan 2 Zgyl[9U ' v] 
2 2 

(74) = 1 + 2 7u7 v 7euev ~ 7u - 7v - 7| u ev 

(1 + 7 U + 7v + 7 Q u©v) 2 

= tan 2 | 

The extreme right-hand side of (JTJ) follows from (J20J) - (EU ■ 

Interestingly, it follows from (|74|) that the Thomas precession angle generated 
by u and ©v, that is, Zgyr[u, Gv], possesses an important hyperbolic geometric 
property. It equals the defect 5 of the gyrotriangle generated by u and v in R^; see 
also [15J pp. 236-237] and the Gyration - Defect Theorem in US] Theorem 8.55, 
p. 317]. 

The gyration gyr[u, 0v] possesses an important gyroalgebraic property as well. 
It gives rise to a second binary operation EB, called Einstein coaddition, given by 
the equation 

(75) u ffl v = u©gyr[u, 9v]v 
which can be dualized into the equation [48j Theorem 2.14] 

(76) u©v = u ffl gyr[u, v]v 

Unlike Einstein addition, which is gyrocommutative, Einstein coaddition is com- 
mutative. Furthermore, it possesses a geometric interpretation as a gyroparallelo- 
gram addition law, and it gives rise to the two mutually dual right cancellation laws 

SB] 

(v©u) B u = v 

(77) 

y J (v ffl u)Qu = v 

Einstein coaddition ffl captures useful analogies with classical results, two of 
which are the right cancellation laws in (|77[) . Another important analogy that Ein- 
stein coaddition captures is associated with the gyromidpoint. Indeed, the midpoint 
M A B shown in Fig. [2] is expressed in terms of the points A and B by the equation 

(78) M AB =A+{-A + B)\ = \{A + B) 

while, in full analogy, the gyromidpoint M A B shown in Fig. [3] is expressed in terms 
of the points A and B by the equation [48] 

(79) M AB = A®{eA®B)®\ = \®{AmB) 
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The right part of Fig. @] raises the question as to whether the composite velocity 
of frame E" relative to frame E is u®v or v©u. The answer is that the composite 
velocity of frame E" relative to frame E is neither u©v nor v©u. Rather, it is 
given by the commutative composite velocity u EH v. Indeed, it is demonstrated 
in [5TJ Chap. 10-Epilogue], and in more details in [351 Chap. 13], that looking 
at the relativistic velocity addition law and its underlying hyperbolic geometry 
through the lens of the cosmological stellar aberration effect leads to a startling 
conclusion: relativistic velocities are gyrovectors that add in the cosmos according 
to the gyroparallelogram addition law of hyperbolic geometry, that is, according to 
the commutative addition uEHv, rather than either Einstein addition u©v or v©u. 



13. Thomas Precession Frequency 

Let us consider a spinning spherical object moving with velocity v of uniform 
magnitude v = ||v|| along a circular path in some inertial frame E. We assume that 
the spin axis lies in the plane containing the circular orbit, as shown in Fig.[7J The 
spinning object acts like a gyroscope, maintaining the direction of its spin axis in 
the transition from one inertial frame into another one, as seen by inertial observers 
moving instantaneously with the accelerated object. Following Taylor and Wheeler, 
we approximate the circular path by a regular polygon of n sides [34] . as shown in 
Fig. [7] for n = 8. In moving once around this orbit the object moves with uniform 
velocity v in straight-line paths interrupted by n sudden changes of direction, each 
through an angle 8 n = 2n/n. 

An observer at rest relative to the laboratory frame E views the motion of the 
object along the polygonal path as the result of successive boosts (A boost being a 
Lorentz transformation without rotation [41]); see Sec. [14j He therefore measures 
a Thomas precession angle e n by which the object spin axis is precessed when the 
object rounds a corner. By (|67|) . this Thomas precession angle e„ is determined by 
the equations 

( 7v -l) 2 sin 2 ^ 



(80) 



( 7 2 + l) + ( 7 3-l)cOsf 



( 7 2 -l) + ( 7v -l) 2 cos^ . 2tt 
sin e n = - — — — — — ~— sin 

(ll + l) + (73-l)cos^ 



By Euler's equation we have 

e le " = cos e n + i sin e r , 



(81) 



where i = -J— I and 



n 



, fi0 x f ( M (7v - I) 2 sin 4> + i{ ill - 1) + (7v ~ I) 2 cos <t>} ■ 
(82) f((p) = — — sm< 



2+( 7 3-l)(l + cos( 



G 
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Figure 7. A regular polygonal path in R 3 as an approximation 
to the Newtonian circular path of a spinning spherical object. The 
change of direction at each vertex of the polygon is 6 n — 2ir/n, 
where n is the number of the polygon sides. Here, n — 8. In the 
limit n —¥ oo, the polygonal path tends to the circular path. A 
spinning spherical object is moving with velocity of uniform mag- 
nitude along the polygonal path. The points A,B,C € M 3 are three 
adjacent vertices of the polygon in the rest (laboratory) frame E. 
When the object moves from A to B it is at rest relative to the 
frame E', and when the object moves from B to C it is at rest rel- 
ative to the frame E". The relationship between the three inertial 
frames E, E' and E" is thus the one shown in Fig. 2] with 9 = n . 
Accordingly, since the object moves in the counterclockwise direc- 
tion, it precesses in the clockwise direction. 

Initially, the spin of the object is vertical when the object moves 
uniformly from A to B. After completing its first closed orbit in 
the counterclockwise direction, the object returns to is original po- 
sition, now moving from A to B with a spin that is precessed in 
the clockwise direction. The initial spin and the final spin for the 
first closed orbit starting at A are shown. 



As the spinning object moves around its polygonal orbit, its spin axis, as observed 
in E, precesses by the Thomas precession angle e„ when it rounds each of the n 
corners of the polygon as shown in Fig. [7] The total angle of precession is thus ne n , 
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represented by the unimodular complex number 

(83) e *"» = (i + /(2E; 

In the limit n — > oo the polygonal path becomes a circular path, and the frame 
of reference in which the center of momentum of the spinning object is momentarily 
at rest is being changed continually. The total Thomas precession is thus the angle 
et given by the equation 

(84) e ltt = lim e m£ " = lim (l + /( — ) 

n— too n— >oc I 77, 

Let g(x), xGl, be the function 

(85) g(x)= lim {l + /(-)}" 

The function f(<f>) is continuous on R, satisfying /(0) = 0. Hence, 

(86) lim /(-) = 

n— >oo n 

for any x £ R. 



Interchanging the limit in (|85|l with a differentiation with respect to x we find 
that the function g(x) satisfies the initial value problem 

g'(x) = f(0)g(x) 

(87) 9(0) - 1 
for x e R. 

The unique solution of the initial value problem (|87[) is 

(88) = e f '(° )x 

Hence, in particular for x — 2tt, it follows from (|M)l . (j85|) and (|88p that 

(89) e l£t = g{2ir) = e 2w/ '(°) 

But, 

(90) /'(0) = 1 

/V 

Hence, by (|89 |) - f|90|) , the Thomas precession angle e t is given by the equation 

7 V -1 



(91) e t = -2tt- 



7v 



The Thomas precession angle et is the angle through which the spin axis precesses 
in one complete circular orbit. It requires, therefore, 2jr/et orbits for the object to 
precess to its original orientation through 2n radians. Hence, if the angular velocity 
of the circular motion of the object is to, then the angular velocity ui t of the Thomas 
precession angle of the object is given by the equation 

(92) oj t = —u = w 

2?r 7 v 
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The quantity ui t in (j9"2")l is the angular velocity of the Thomas precession angle et 
of a particle that moves in a circular orbit with angular velocity u>. Eq. (|92|) relates 
the angular velocity uj t of the Thomas precession angle e t to its generating angular 
velocity us. It demonstrates that the angular velocities iv t and u are oppositely 
directed, as shown graphically in Fig. 01 

If the magnitude of the velocity v and the acceleration a of the spinning object 
are v and a then its angular velocity is given by the equation u = a/v. Hence, the 
angular velocity uit of the Thomas precession angle et is given by the equation 

(93) u t = 

7 V v 

Taking into account the direction of the Thomas precession axis and the velocity 
and the acceleration of the spinning object, and noting (|10[) . the Thomas precession 
angular velocity u> t in (|93|) can be written as a vector equation, 

(94) Wt = 7^1axv = _T^axv 

7v v 1 + 7v c 

The coordinate axes in the rest frame of any body in torque-free, accelerated 
motion precesses with respect to the laboratory axes with an angular velocity uj t is 
given by ([94]). Since 7 V /(1 + 7 V ) = 1/2 + (l/8)(u 2 /c 2 ) + . . ., the angular velocity 
uit of the resulting Thomas precession, for the case when v = ||v|| « c, is given 
approximately by the equation 

. . 1 a x v 

(95) uh = ~ — 

According to Herbert Goldstein [T5] p. 288], u>t in (|93|) is known as the Thomas 
precession frequency. 

The Thomas precession frequency (|9"5")) involves the famous factor 1/2, known as 
Thomas half. The experimental significance of this factor is well known: The spin- 
ning electron of the Goudsmit-Uhlenbeck model gives twice the observed precession 
effect, which is reduced to the observed one by means of the Thomas half [57] . 



14. Thomas Precession and Boost Composition 

Einstein addition underlies the Lorentz transformation group of special relativity. 
A Lorentz transformation is a linear transformation of spacetime coordinates that 
fixes the spacetime origin. A Lorentz boost, -B(v), is a Lorentz transformation 
without rotation, parametrized by a velocity parameter v g R^. The velocity 
parameter is given by its components, v = (vi, 112,1)3), with respect to a given 
Cartesian coordinate system of M.^.. Being linear, the Lorentz boost has a matrix 
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representation, which turns out to be |23j . 



S(v) 



(96) 



( 7v 


c 2 7v w i 


C 2 7v w 2 


C 2 7v«3 \ 


7v u i 




— 2 T 2 


— 2 "Y 2 


7v«2 


— 2 X 2 , 

C ^1^2 


1 + c- 2 7 l lW 2 


— 2 7^ 


\7v^3 


— 2 7-, 2 , 


—2 7^ 
C 7v + x «2W3 


l + c- 2 ^T«3V 



Employing the matrix representation (|96|) of the Lorentz transformation boost, 
the Lorentz boost application to spacetime coordinates takes the form 



(97) 



fl(v) 



S(v) 















X 2 






w 




U/ 



where v = («i,«2,U3)* € R 3 , x = ( x ii ^2, £3)* € R 3 , x' = (x' x , x' 2l x' 3 Y £ R 3 , and 
t, t' € R, where exponent t denotes transposition. 

A 2-dimensional boost is obtained in the special case when v$ — X3 = in (|96p 
(|97|l. For simplicity, the boosts of inertial frames in Fig. 0] are two dimensional 
boosts, and time coordinates are not shown. In this figure, the spacetime coordinate 
systems E, £' and E" (only two space coordinates are shown) are related by boosts. 
Specifically, in Fig. [H 

(1) the application of the boost B(u) to the spacetime coordinate system E 
gives the spacetime coordinate system £', 

(2) the application of the boost B(v) to the spacetime coordinate system E' 
gives the spacetime coordinate system E", 

(3) and the application of the boost B(u®v), or £?(v©u), to the spacetime co- 
ordinate system E, preceded, or followed respectively, by a Thomas preces- 
sion (see (|101j) - (|102p in Theorem [7] below) gives the spacetime coordinate 
system E". 

The Lorentz boost (l9"6"|) - (|9"T|) can be written vectorially in the form 



/A / 7u(i+^u-x) \ 

(98) B(u) U = , 
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Rewriting with x = vt e M. 3 , u, v e K 3 C M 3 , we have 



B(n)[l )=B(u) , , , 



(99) 



7- 



Tv 



(u®v)t , 



The equations in (J99| reveal explicitly the way Einstein velocity addition under- 
lies the Lorentz boost. The third equality in (|99| follows from (|9a| and (|2]l. 

In general, the composition of two boosts is equivalent to a single boost preceded, 
or followed, by the space rotation that Thomas precession generates, as we see from 
the following theorem: 

Theorem 7. (The Boost Composition Theorem): Let u,v,w 6 K 3 be rela- 
tivistically admissible velocities, let x = wt, t > 0, and let B(u) and -B(v) be two 
boosts. Furthermore, let Gyr[u, v] be the spacetime gyration of space coordinates, 
given by 



( 10 °) G ^[ U ' V ](xiwt) : -((gyr[u!v]w)t 



Then, boost composition is given by each of the two equations 
(101) B(u)B(y) = J B(uffiv)Gyr[u,v] 



(102) B{u)B{v) = Gyr[v,u]£(v®u) 



Proof. We will show that (|101l) follows from the gyroassociative law of Einstein 
addition and that that (jl02[) follows from (|101[) and the gyrocommutative law of 
Einstein addition. 



:S2 



On the one hand we have the chain of equations below, which are numbered for 
subsequent derivation: 



L'(ui/J(v)( B(u)B(v)(^ 



(2) / Tvewj 



S(u) 



7. 

v 



(v©w)i 



= B(u) ' 



(103) 



(v©w)*' 



(4) / 7 uffi(vffiw ) ^/ 

7u S (v S w) { u@ ( v@w )W 



(5) / 7 uffi(vffi w) ■ 

7u y ) {ti9(v$w)}t / 



(6) / 7 (uffiv)ffiByr [ u , v ] w ^ 

7 (u ffi v) e gyr [ u . v ]w {( u@v ) @gyr [ U) y] W }t 



Derivation of the numbered equalities in fjl03|) follows: 

(1) Follows from the definition x = wi. 

(2) Follows from Q] by a boost application to spacetime coordinates according 

to HMD- 

(3) Follows from [2] by the obvious definition 



(104) t' = 



(4) Follows from [3] by a boost application to spacetime coordinates according 

to dag) . 

(5) Follows from [4] by the substitution of (|104l) for t'. 

(6) Follows from [5] by the gyroassociative law of Einstein addition. 
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On the other hand we have the chain of equations below, which are numbered 
for subsequent derivation: 



J B(u©v)Gyr[u,v] 



(105) 



(2) 



(3) 



(4) 



£(u©v) 
S(u©v) 



gyr[u, v]x 
(gyr[u,v]w)t 

T(u8v)egyr[u,v] 



-t 



Tgyr[u, V ]v 

Vwter[u,v)w {(u0v)©gyr[u,v]w}t 

\ Tgyr[u,vlw 



7 (uev)egyr[u.v]w £ 

7 (uev)egyr[u,y] 

7 W 



{(u©v)©gyr[u, v]w}i 



Derivation of the numbered equalities in (|105l) follows: 



(1) Follows from the definition of the spacetime gyration Gyr[u, v] in terms of 
the space gyration gyr[u,v] in (llOOp . 

(2) Follows from [1] by definition, x = wt. 

(3) Follows from [5] by a boost application to spacetime coordinates according 
to (UlI)- 

(4) Follows from [3] by the identity 7 W = 7 gyr r uv i w that, in turn, follows from 
the definition of gamma factors in ([3]) along with the invariance (|33l) of 
relativistically admissible velocities under gyrations. 



The extreme right-hand sides of the chain of equations (jl03[) and (|105[) are iden- 
tically equal. Hence, the extreme left-hand sides of (| 103[) and (|105[) are equal for 
all spacetime events (i,x)*, t € E, x = wt, w £ Rj?, thus verifying (|101[) . 
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In order to verify (I102p . let us now consider the chain of equations below, which 
are numbered for subsequent derivation: 

Gyr[v,u]B(veu) Gyr[v, u]£(v®u) 



(2) 

= Gyr[v,u] 



T(veu)ew * 
-{(v©u)©w}i 



7w 

,v u, ®"-gyr[v,u]{(v©u)©w}t 



-t 

7 W 



(106) 



(4) / 7gy,[u,v]{(veu)il)w) ^ \ 
7 gyr[ U , V ]w 

Vr[u,v| (( .e uM w) gyr | V| u ]{(vffiu)ffiw}i 

'gyr[u,v]w / 

(5) / 7 (uev)e8yr[u.y]w ^ 

T E yr[ U , V ]w 

7 (uffi v)e g yr[u,v]w | ( U@V ) @gyr [ Ui V ] 



7, 



gyr[u, V ]v 



(6) 



= B(uffiv) 

(7) 



gyr[u, v]wf 



B(uffiv) ( r , 
y ^gyr[u,v]x 

X ft 
= = B(uffiv)Gyr[u,v] ( 

The chain of equations (j 106[) is valid for all spacetime events (t, x)*, t € E, x = wi, 
welj, thus verifying (fT02|) . 

Derivation of the numbered equalities in (|106l) follows: 

(1) Follows by definition, x = wt. 

(2) Follows from Q] by a boost application to spacetime coordinates according 
to ([99]). 

(3) Follows from[2]by the definition of the spacetime gyration Gyr[v, u] in terms 
of the space gyration gyr[v, u] in (|100|) . 

(4) Follows from[3]by the identity 7 W = 7 gyr [ u , v ] w > u, v,w e that, in turn, 
follows from the definition of gamma factors in ([3]) along with the invariance 
(|33[) of relativistically admissible velocities under gyrations. 

(5) Follows from|4]by the linearity of gyrations along with the gyrocommutative 
law of Einstein addition. 
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(6) Follows from [5] by a boost application to spacetime coordinates according 

to (35). 

(7) Follows from [5] by definition, x = wt. 

(8) Follows from[7]by the definition of the spacetime gyration Gyr[v, u] in terms 
of the space gyration gyr[v, u] in (jlOOp . 



The Boost Composition Theorem [7] and its proof establish the following two 
results: 

(1) The composite velocity of frame S" relative to frame E in Fig. |4] may, 
paradoxically, be both u©v and v©u. Indeed, it is u©v in the sense that 
E" is obtained from E by a boost of velocity u©v preceded by the gyration 
Gyr[u, v] or, equivalently, it is obtained from E by a boost of velocity v©u 
followed by the gyration Gyr[v, u]. 

(2) The relationships (|101[) - (|102[) between boosts and Thomas precession are 
equivalent to the gyroassociative law and the gyrocommutative law of Ein- 
stein velocity addition as we see from the proof of Theorem [7] 

In view of these two results of the Boost Composition Theorem, the validity of the 
Thomas precession frequency, as shown graphically in Fig. |4j and the relationship 
between the Thomas precession angle e and its generating angle 9 stem from the 
gyroassociative law of Einstein velocity addition. Hence, in particular, the result 
that e and 9 have opposite signs is embedded in the gyroassociative law of Einstein 
addition. In the next section we will present a convincing numerical demonstration 
that interested readers may perform to determine that, indeed, e and 9 in Fig. |4] 
have opposite signs. 

15. Thomas Precession Angle and Generating Angle have Opposite 



As in Fig. [4j let e and 9 be the Thomas Precession Angle and its generating 
angle, respectively. As verified analytically, and as shown graphically in Fig. [4j the 
angles e and 9 are related by (|66[) and, hence, they have opposite signs. 

Without loss of generality, as in Fig. |4j we limit our considerations to two space 
dimensions. Let u, v € M.^ be two nonzero relativistically admissible velocities with 
angle 9 between their directions, as shown in Fig. 2] Then, they are related by the 
equation 



Let w e Mj be the velocity of an object relative to frame E" in Fig. |U Then, 
the velocity of the object relative to frame E in Fig. @] is 



□ 



Signs 



(107) 




(108) 



u©(v©w) = (uffiv)ffigyr[u, v]w 
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so that the velocity w of the object is rotated relative to E by the Thomas precession 
gyr[u, v], which corresponds to the rotation angle e given by (|66|) . Hence, 

, , r i /cose -sine\ 

(109) gyr u, v w = . w 
v ' Ume cosey 

where e is given by (|66|) . 

Substituting v from (|107[) into (I109[) . we obtain the equation 

,„.,_.. r ll v ll /cos 6* — sin#\ . /cose — sine\ 

(110) gyr u, -— - . a fl u w = . w 
' L u \sm6> costfy J Vsme cosey 

In (jllOp . 9 is the angle shown in the left part of Fig. [H which generates the 
Thomas precession angle e shown in the right part of Fig. [4j where e is determined 
by 9 according to and, hence, where 9 and e have opposite signs. The validity 
of (|110l) can readily be corroborated numerically. The numerical corroboration of 
the validity of (I110[) . in turn, provides a simple way to convincingly confirm our 
claim that indeed 9 and e have opposite signs. 
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